Exact Parent Hamiltonian for the Quantum Hall States in a Lattice 
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We study lattice models of charged particles in uniform magnetic fields. We show how longer 
range hopping can be engineered to produce a massively degenerate manifold of single-particle 
ground states with wavefunctions identical to those making up the lowest Landau level of continuum 
electrons in a magnetic field. We find that in the presence of local interactions, and at the appropriate 
filling factors, Laughlin's fractional quantum Hall wavefunction is an exact many-body ground state 
of our lattice model. The hopping matrix elements in our model fall off as a Gaussian, and when 
the flux per plaquette is small compared to the fundamental flux quantum one only needs to include 
nearest and next nearest neighbor hoppings. We suggest how to realize this model using atoms 
in optical lattices, and describe observable consequences of the resulting fractional quantum Hall 
physics. 
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The interplay between periodic potentials and mag- 
netic fields is an important topic |l|-|5| . In the tight bind- 
ing limit, the lattice broadens the Landau levels into a 
series of finite bandwidth "Hofstadter bands" which can 
be represented as a self-similar fractal. Since the original 
band-gaps persist, the integer quantum Hall effects are 
robust against the lattice. The split degeneracy, how- 
ever, invalidates many of the analytic arguments used 
to explain the fractional quantum Hall effect [a-ISI], and 
questions remain about the nature of the interacting sys- 
tem. Here, by adding longer range hoppings to a Hub- 
bard model, we produce a Hamiltonian for which sev- 
eral Hofstadter bands coalesce into a single degenerate 
manifold. Adding local repulsion between the particles, 
we show that at appropriate filling factors the Laughlin 
wavefunction becomes an exact ground-state. 

In a uniform magnetic field, the most general hopping 
Hamiltonian on a two-dimensional square lattice is 



H = y^J(zj,Zk)a\ak, 



(1) 
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J(z 3 ,z k ) = W{z)e^ ,2) ^ z '- z 'o z )*, 



where the position of the j'th lattice site is written in 
complex notation as Zj = Xj + iyj , and z = Zk — Zj . The 
operators Oj annihilate an atom at site j. The phase 
factor (zjz* — z*z) tfi = 2i (xj-y — yjx) <fi, corresponds to 
a uniform magnetic field in the symmetric gauge, with 
flux <f> through each plaquette. This flux is only defined 
modulo 1, and having a full flux quantum through each 
plaquette is gauge equivalent to no flux. We will explic- 
itly assume < <j> < 1, and take <f> = p/q to be the ratio of 
two relatively prime integers. If one chooses W to be —t 
for nearest neighbors and zero otherwise, one reproduces 
the Hofstadter spectrum [l| . We show that if instead we 
choose 



^)=(xG(z)e^[ (1 "*l ! 
G (z) = (-lf +y+xy , 



(2) 



the lowest p Hofstadter bands collapses to a single fully 
degenerate Landau level. Although we work in the sym- 
metric gauge A = (B/2) (xy — yx), converting our re- 
sults to other gauges is trivial: under a gauge transfor- 
mation A(r) — > A(r) + VA(r) and Cj — > Cje lA( - r ^. The 
flux is measured in units of 0o == h/e, where h is Planck's 
constant, and e is the electric charge. Our derivation of 
this Hamiltonian is similar to one used by Laughlin [10( 



and subsequently corrected/extended in [ll|, |12|. The 
paradigm of creating a parent Hamiltonian for which a 
desired quantum state is an exact eigenstate has been 
fruitful in a number of other spin models 13Wl5J|. We 
work in units where t = 1 . A similar construction can be 
defined for triangular lattices |16j |. 

Our results have deep implications. First and fore- 
most, they provide an exact equivalence between the 
lowest Landau level in the continuum and in a realis- 
tic lattice system. This equivalence is unexpected, and 
can be further exploited. For example, it provides an 
avenue for robust lattice calculations of continuum quan- 
tum Hall systems, and investigations of lattice fractional 
quantum Hall states [4|,ll7|. Similarly it provides a means 
for lattice experiments to emulate an important contin- 
uum problem (18J . 

Massive degeneracies, such as the ones we found here, 
are related to symmetries. Further theoretical work on 
this model may reveal these symmetries. Our model also 
provides a convenient inroad towards a semiclassical un- 
derstanding of the Hofstadter spectrum, and connections 
to magnetic breakdown in real materials |19j |. 

The most promising experimental realization of our 
model is in optical lattices [18(. Optical lattice exper- 
iments can study both bosonic and fermionic quantum 
Hall states, and allow us in principle to study much 
larger fluxes than can be achieved with real magnetic 
fields. The gauge potential in the optical lattice system 
can be created in a number of ways: time-varying hop- 
ping elements |2|, lattices with multiple sets of minima 



1201 , coherent Raman scattering [2l| and rotation |5|, |22| - 
l26j . Further, optical lattice systems allow us to directly 
tune the hopping amplitudes between nearby sites. Long 
range hopping is difficult to arrange, but in our model 
J falls off as a Gaussian, and in the limit of small </> it 
suffices to include only nearest and next-nearest neighbor 
hopping. The ratio of these hopping matrix elements can 
be controlled in an experiment by adjusting the shapes 
of the barriers between those sites. One practical scheme 
would involve adding an additional array of shallow wells 
displaced by half a lattice spacing in both the x and y 
directions. Integrating out these shallow sites will renor- 
malize the nearest and next-nearest neighbor hoppings. 
A second scheme would be to divide the original lattice 
into two sublattices. Separating the sublattices in the 
z-dircction will attenuate the nearest neighbor tunelling 
while leaving the next-nearest neighbor matrix element 
largely unchanged. Figure [T] compares the energies of 
the single particle eigenstates of Eq. (P), using the W 
in Eq. ([2]), as well as truncating to only nearest neigh- 
bors or next-nearest neighbors. As one can see, even for 
<f> = 1/3, the next-nearest neighbor hopping already re- 
duces the bandwidth to O.lt. 

Another possible realization of Eq. ([2]) would be 2D 
electron gases in a superlattice 27H30J. Hopping ampli- 
tudes can be tuned through altering the device structure. 
These systems also naturally include long-range Coulomb 
interactions, which are absent in trapped neutral atoms, 
and lead to richer many-body physics. 

Not only does this Hamiltonian produce a macroscopi- 
cally degenerate manifold of single particle ground states, 
but this manifold is spanned by wavefunctions of the form 



i-„ ( :,) : ;./>•„: : : zf exp ( — — | Zj 1 2 



(3) 



all with energy e = — I. Remarkably, this is the same 
structure as the continuum problem, where the LLL is 
characterized by the same degenerate set of single particle 
states. To prove this result, we write 



{j\H\<M 

(Mn) 



= E G w 



(zj + z) r 



\z\ — 7T0Z*Z 



2^0 



We then appeal to the singlet sum rule [lOj, |31| 



fc(c)^E eCZG ( 



z)e 2 



M 2 



Vc, 



(4) 



(5) 



where the sum is over all z = n + im with integer n and 
m. By taking any number of derivatives with respect to 
c one finds 



£/(*)G(z)e- 



-N ; 



0, 



(6) 



for any entire function / (z) that diverges sufficiently 
slowly as \z\ — > oo 32j. Since we do not include the 
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FIG. 1: All single particle eigenvalues for the hopping Hamil- 
tonian in Eq. JT} with <f> — 1/3 on a 12 x 12 lattice with 
periodic boundary conditions. The index n labels eigenvalues 
from smallest to largest. The white disks use hopping ma- 
trix elements given by Eq. J5J, the black disks are the same 
model with only nearest and next nearest neighbor hopping, 
and the grey boxes have only nearest neighbor hopping (the 
Hofstadter Hamiltonian). Energies are all measured in units 
of t. The energies of the white and black disks are nearly in- 
distinguishable, and the white disks obscure some of the black 
ones. The lowest 1/3 of the white disks are all degenerate. 



z = term in Eq. (|4]), one immediately finds that the 
right hand side is simply -1, proving that the LLL wave- 
functions ([3]) are degenerate eigenstates. No analogous 
argument works for the higher Landau level wavefunc- 
tions, which involve powers of both z* and z. 

Given that the wavefunctions in (j3|) are identical to 
those of the continuum problem, the total number of de- 
generate states per unit area is the same as in the contin- 
uum; this results in <pN s LLL wavefunctions in a region 
containing N s lattice sites. Thus <f> is the fraction of all 
single particle states which reside in the LLL. Taking 
4> = p/q> the standard Hofstadter problem yeilds q dis- 
tinct bands. Thus, as we confirm numerically, our LLL 
must be made from the lowest p of these. This p-fold 
collapse is consistent with the relationship between the 
Chern numbers of the Hofstadter bands, and that of the 
LLL [33J]. 

For </) > 1/2 it is natural to also consider the Hamilto- 
nian formed if one replaces 4> in Eq. (2) with 1 — <f> and 
leaves equation (1) unchanged. Due to the periodicity 
in <f) of lattice models, this gives a Hamiltonian with the 
same absolute flux per plaquette, however it is clearly a 
distinct Hamiltonian, with shorter range hopping. This 
alternative Hamiltonian yields states analogous to ((3|), 
but with z replaced by z* and a degeneracy o£(l — (j))N s , 

The massive ground state degeneracy of our system can 
be lifted by interactions. Since our model reproduces the 
continuum lowest Landau level, we can simply use those 
results [3, y, |34{ . On-site repulsion in the lattice is equiv- 
alent to point interactions in the continuum. Consider 
for example the interacting Hamiltonian 



H = 



3^k,<T 



J{zj,zu)a] (T aka + 



E 



u„ 



- a ja a ja' a j<T' a J<r( 7 ) 



where U aa i > is the on-site interaction energy of par- 



tides with spin states a and a' . Any LLL wavefunction 
which vanishes when two particles coincide is a ground 
state of this Hamiltonian. Due to the structure of the 
LLL, there is a maximal atomic density for which this 
occurs. For single component bosons the highest density 
ground state is the v = 1/2 Laughlin state, 



JY 



iV 



i/>(zi,--- ,z N ) = Y[ ( Z i ~ z if II ' 



(8) 



K] = l 



For 2-component fermions it is the ferromagnetic "111" 
state 34[ . At fixed density, these states are unique up to 
topological degeneracies. The other Laughlin states and 
more exotic quantum Hall states are also ground states if 
v < 1/2, however they are not unique if all interactions 
are local. 

Longer ranged interactions (as found in electronic sys- 
tems or in optical lattice experiments with dipolar gases 
[26|) will typically lift the degeneracy entirely |35|. The 
subsequent analysis can be quite involved, but most con- 
tinuum arguments will carry over. 

Our construction is readily extended to a finite system 
with magneto-periodic boundary conditions, 



i/j (z + nL + imL) = ijj (z) e 



iTT(f>L(ny—mx) 



(9) 



There one replaces the polynomials in Eq. (|3]) with ap- 
propriate products of Gaussians and Jacobi theta func- 
tions 36j, |37|. One also replaces J(zj,z) in Eq. (JTJ) by its 
magneto-periodic extension 

J L (Z 3 ,Z) = Y^ J ( Z 3' Z + R ) ex P (o ( Z 3 R * ~ Z 3 R )<P) > 
R 



where the sum is over all R = nL + imL for integer 
n.m. This finite system is amenable to numerical cal- 
culations. To invoke the singlet sum rule in a periodic 
geometry, one must simply merge the sums on z and R 
into a single sum over all z =/= 0. The phase factors from 
the magnetoperiodicity of ip (z) and J cancel each other. 
In the finite system with magneto-periodic boundary 
conditions the v = 1/p Laughlin state is 

M M 

k<j j = l 

M> cm = n^(^(z-zo). (11) 

The center of mass coordinate is Z = ^\. Zj, and Xjk = 

0i(tt(*j -z k )/L), with 9i(z) = E„(-l)"" 1/2 e-(" +1 /2) 2 
e jz(2n+i)_ There are p parameters Zi which represent the 
location of the center of mass zeros. In the continuum 
system there is a symmetry which causes the energy to 
be independent of how these are chosen. The space of 
degenerate states is spanned by p orthogonal wavefunc- 
tions. In most lattice models this symmetry is broken, 
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FIG. 2: First 100 eigenvalues for 4 particles on a 4 x 4 lattice 
with periodic boundary conditions, <jf> = 1/2, and hard core 
repulsion. The two states at e = —4 are Laughlin states (111ft : 
the degeneracy stems from the toroidal geometry. There is a 
distinct energy gap of 0.566t n n to the lowest excited states, 
where t nn is the nearest neighbor hopping amplitude. 



and the degeneracy is lifted. Since Eq. (fTTj) is made up of 
lowest Landau wavefunctions, in our model the degener- 
acy persists. In Fig. [2] we confirm this degeneracy via an 
exact diagonalization calculation for 4 bosons on a 4 x 4 
lattice with p = 2 and hard-core repulsion. 

Our results gives some insight into recent calculations 
of Sorensen et al. [2J, [26|. They investigated the stan- 
dard Bose- Hubbard model with nearest neighbor hopping 
and a uniform magnetic field. Fixing the filling factor at 
v = 1/2, they found that when (j> became of order 0.2 the 
overlap between the exact ground state and the p = 2 
Laughlin state (fTTj) begins to rapidly decrease. The char- 
acteristic range of hoppings in our model increases with 
(10) (f> - and near <j) — 0.2 the next nearest neighbor matrix 
element starts to become significant. 

Since we advocate using cold atoms to investigate this 
physics, it is important to understand how FQH physics 
would manifest itself in those systems. Although most 
difficult, the most exciting observations would be ones 
which probe the braiding properties of the excitations 
[38l l39j . These states also have definite signatures in 
Bragg spectroscopy [3j. The most robust probe, how- 
ever, is an analog of the vanishing longitudinal resistance 
seen in solid state systems - namely the incompressibil- 
ity of the fractional quantum Hall states [a, |40j . This 
incompressibility is readily observed in trapped systems, 
where the chemical potential (and hence the filling fac- 
tor) varies slowly in space. As is caricatured in Fig. [31 
the equation of state n(/x) has a series of plateaus corre- 
sponding to the filling factor taking on integer fractions. 
Within the local density (Thomas- Fermi) approximation, 
the density pro file of t he trapped cloud will display these 
same plateaus 



-431]. The width of these plateaus is 
set by the gap to single particle excitations in the frac- 
tional quantum Hall states. As shown in figure [51 in the 
hard core limit the gap in a 4 x 4 lattice at v = 1/2 and 
<p = 1/2 is 0.566 t nn . This should be compared to the 
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FIG. 3: Schematic plot of (n) vs. p, for lattice bosons de- 
scribed by Eqs. (|1I2[I . with hard core interactions added. The 
steps correspond to incompressible particle/hole fractional 
quantum Hall states at v — 1/2. This structure will be visible 
in the density profile of a trapped gas. Similar structure will 
be seen with Fermions, but with plateaus at fillings with odd 
denominators. 



bandwidth V ~ 4t nn . As [i goes from to V the density 
goes from zero to one. One therefore expects that the 
v = \j1 plateau will occupy roughly 1/8 of the cloud. 
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